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Manipulation of the effects of consciousness by external influence on the 
human brain is considered in the context of the nonlinear dynamical modeling 
of interaction between automatic and conscious processes. 



In previous papers [1,2] an approach to nonlinear dynamical modeling of 
interaction between automatic (A) and conscious (C) processes in the brain was 
presented. The idea is to use quaternion field with real and imaginary compo- 
nents representing A - and C - processes. The subjective C - experiences were 
divided into three major groups: sensations (S), emotions (E) and reflections 
(R). Note, that subjective S should be distinguished from the automatic sen- 
sory input into the neuron system of the brain. The A - C interaction is due to 
the nonlincarity of the system. This approach was illustrated on the nonlinear 
equation for the current density in the cortex. The nonlincarity is determined 
by the sigmoidal firing rate of neurons. Perspective for testing of this approach 
were also indicated as well as some more general approaches [1,2]. 

For the purpose of medical and other possible applications it is interesting 
to include an external electromagnetic (EM) influence in this modeling. In a 
laboratory setting a specially equipped helmet can produce designed nonhomo- 
geneous or homogeneous excitations in the brain. On another hand, suppose 
we want to pacify a group of terrorists (!) by using a strong EM radiation 
with the wavelength much larger than the size of their brains. In this case the 
excitation will be approximately homogeneous. We start with the homogeneous 
case which is more simple mathematically and gives some insight into general 
situation. 

The model equation for the average (spatially uniform) current density a(t) 
perpendicular to the cortical surface has the form [1,2]: 
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Here k is the relaxation coefficient, a(t) is the average sensory input, / represents 
the sigmoidal firing rate of neurons [for example, f(a) — tanh(a)], components 
4> p represent the indicated above (S, E, R) - effects and summation is assumed 
on repeated subscripts from 1 to 3. The quaternion imaginary units i p satisfy 
conditions: 
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where s pqr is the unit antisymmetric tensor and 5 pq is the unit tensor. Formula 
(2) is a compact form of conditions: i\ = i\ = i% = — 1, i\i2 = — «2*i = h, 
= — «3«2 = ii, hh = —iih = h- Equation (1) is obtained by using the 
quaternion q = a + i p ip p instead of a in order to describe the A - C interaction. 
The additional term ip in (1) represents the external EM excitation. Equation 
(1) is the real part of the equation for the quaternion [1,2]: 

^+kq = f(q + cr) + <P ((la)) 
For il) p from (la) we have equations: 

-jg- + fcV P = Im p {/(a + <j + i q i> q )}, p = 1, 2, 3 ((3)) 

where Im p {/} = — Ke{fi p }. Note, that so-called extra-sensory effects (if they 
exist) can be included in this approach by assuming that a is a quaternion: 
a ==>- a + i p s p , this will produce shift ip p => ip p + s p in the nonlinear terms in 
(1), (3) and below in (5), (6). 

Let us consider typical f(a) — tanh(a). Simple algebra gives [2]: 

, . . sinh(2a) + 7 sin(2-0) ,2 ,2 , , -1 2 , „,,, 

tanhW = cosh(2 Q)+ cos ( 2Jr * " ^ h ^ > 3 1 m 
Using (4) with shift a a + <j, we rewrite (1) and (3) explicitly: 

da I ka _ sinh f 2 ( a + g )] v6 ((5)) 

8t +ka ~ cosh[2(a + a)} + cos(2^) + U >> 

&I>p . t 1 V'pV'~ lsm (2^) 

— ii + fcV = rr- ——,p= 1,2,3 (6 

ot p cosh[2(a + er)J + cos(2?/>) 

Some general conclusions can be made without solving these equations. 
Firstly, if ip p {0) = than ip p (t) = (unless s p ^ 0). Secondly if ^ p (0) ^ 0, 
than evolution i/> p (i) can be manipulated by using sensory input a(t) and EM 
excitation <j>(t). Thirdly, the nonlinearity of the system suggests that the effi- 
ciency of such manipulation depends not only on the amplitudes of <j(t) and 
(j)(t) but also on the shape of these functions (spectral content). 

For the case of spatially nonuniform a(t, x), ip p (t, x), a(t, x) and x) we 
can use more general equations, which include typical propagation velocity of 
signals in the neuron system of the cortex v. Time differentiation of (la), simple 
algebra and addition a term with the two-dimensional spatial Laplacian A gives 
[1,2]: 



_| + {k + m )_| + {km _ v 2 A)q = (m + -)f(q + °) + f t ((7)) 

where m is an arbitrary parameter (see below). Real and imaginary projections 
of (7) give equations for a and tp , which are generalizations of (1) and (3). If we 
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put ip p = and <f> = 0, than equation for a will be similar in spirit to equations 
used for interpretation of EEG and MEG spatial patterns ( see recent paper 
[3] and references therein). In this context we have parameters: k ~ m ~ v/l, 
where / is the connectivity scale. For f(a) = tanh(a) the nonlinear term f(q+cr) 
in (7) has the same projections as in (5) and (6). 

The obtained in this letter equations can be used for numerical experiments 
and for comparison with corresponding laboratory experiments. 
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